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Abstract 

Bethe Ansatz solutions of the open spin-^ integrable XXZ quantum spin chain at 
roots of unity with nondiagonal boundary terms containing two free boundary param¬ 
eters have recently been proposed. We use these solutions to compute the boundary 
energy (surface energy) in the thermodynamic limit. 


In memory of Daniel Arnaudon. 


1 Introduction 


While the solution of the open spin-i XXZ quantum spin chain with diagonal boundary 
terms has long been known [ 111210 ], the solution of the general integrable case, with the 
Hamiltonian mini 


H = Ho + - sinh y coth a_ tanh (d-al + cosech a_ sech f3- (cosh O-a^ + i sinh O^af) 


— coth Q!+tanh/9+(T^ + cosecha+sech/?+( cosh 0 + 0 '^ + z sinh 0 + 0 '^) , (1.1) 


where 



( 1 . 2 ) 


n=l 


(which contains also nondiagonal houndaxj terms) has remained elusive. Here ^ ^ire 

the standard Pauli matrices, y is the bulk anisotropy parameter, a±,j3±, 9± are arbitrary 
boundary parameters, and N is the number of spins. 

Progress has recently been made on this problem. Indeed, a Bethe Ansatz solution is 
now known piziEiin] if the boundary parameters obey the constraints 

1 — €-2 

a- + ei(3- + 620 + + 63 / 3 + = eo(6'_ - 9+) +yk^ -^—ztt mod (2z7r), 616263 = +1, (1.3) 

where 6* = ±1, and k is an integer such that \k\ < — 1 and — 1 + fc is even. Finite-size 

effects for this model and for the boundary sine-Gordon model |3] have been computed on 
the basis of this solution [1211111 . (Related results have been obtained by different methods 
in IEI-) Many interesting further applications and generalizations of this solution have also 
been found (see, e.g., 0). 

Additional Bethe Ansatz solutions with up to two free boundary parameters have been 
proposed in [HI eg. Completeness of these new solutions is straightforward, in contrast to 
the case ca i. A noteworthy feature of the solution [121 is the appearance of a generalized 
T — Q relation of the form 


t{u) Qi{u) = Q2iu') + Q2 {u'') , 
t{u) Q2 {u) = Qi{uf + Qi{u”), 


(1.4) 
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involving two Q-operators, instead of the usual one M- However, unlike the case (m, 
these new solutions hold only at roots of unity, i.e., for bulk anisotropy values 


V 


in 

p + 1 ’ 


(1.5) 


where p is a positive integer. 

The aim of this paper is to use the new solutions [HIISl to investigate the ground state 
in the thermodynamic {N —> cx)) limit. For dehniteness, we focus on two particular cases: 


Case I: The bulk anisotropy parameter has values (EH) with p even; 
the boundary parameters P± are arbitrary, and a± = rj, 9^ = 0 m ( 1 - 6 ) 


Case II: The bulk anisotropy parameter has values (USD with p odd; 

the boundary parameters a± are arbitrary, and (3± = = 0 na (1.7) 


We also henceforth restrict to even values of N. For each of these cases, we determine the 
density of Bethe roots describing the ground state in the thermodynamic limit, for suitable 
values of the boundary parameters; and we compute the corresponding boundary (surface) 
energies. ^ We find that the results coincide with the boundary energy computed in jTU] for 
the case dn, namely, ^ 


'^boundary 


Sin /X / 

2p J-c 


duj 


sinh((i/— 2)a;/4) 1 

2cosh(a;/2)l 2 sinh(z/a;/4) 2 


+ (2 ^^sinh((z/— | 2 a± — l|)a;/ 2 ) ^ sinh(a;/ 2 ) cos( 6 ±a;) ] 


sinh(i/a;/ 2 ) 


sinh(i/a;/ 2 ) / 


1 • / 5 1 

+ - smp cot(/ia±)--cos/i , 


( 1 . 8 ) 


where 


TT 


TT 


/i = —IT] = 


P + 1 ’ 


U = — = p+ I , 
/i 


a± = i/j,a± , f3± = jnb± , 


(1.9) 


and sgn(n) = -j^ for n 7 ^ 0 . 

The outline of this paper is as follows. In Section |21 we briefly review the transfer matrix 
and its relation to the Hamiltonian (dD. In Sections El and m we treat Cases I (USD and 11 
dd, respectively. We conclude in Section El with a brief discussion of our results. 

^For the case of diagonal boundary terms, the boundary energy was first computed numerically in 1^, 
and then analytically in m 

^Here we correct the misprint in Eq. (2.29) of Eli already noted in m 
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2 Transfer matrix 


The transfer matrix t{u) of the model is given by |S] 

t{u) = tiQ K^{u) To{u) Kq{u) fo{u), 

where To(m) and To(m) are the monodromy matrices 

To (m) = Ron (m) ■ ■ ■ Tqi (m) , Tq (n) = Tqi (n) ■ ■ ■ Ron (u) , 
and tro denotes trace over the “anxiliary space” 0. The R matrix is given by 

\ 

n QinVi 7/ ciinln n 0 

R{u) = 


( 2 . 1 ) 


( 2 . 2 ) 


f sinh(M + 7]) 0 0 0 

0 sinh u sinh rj 0 

0 sinh T] sinh u 0 

\ 0 0 0 sinh(M + rj) J 


(2.3) 


where r] is the bnlk anisotropy parameter; and K^{u) are 2x2 matrices whose components 
are given by 

= 2 (sinh a_ cosh/d_ coshM + cosh a_ sinh/9_ sinh m) 

^22 ('^) = 2 (sinh a- cosh f3- cosh u — cosh a_ sinh /?_ sinh u) 

— e^"sinh2n, (n) = sinh 2 m , (2.4) 


and 


K+{u) = K {-u- , (2.5) 

where , jSzp , 6*:^ are the boundary parameters. 

For M = 0, the transfer matrix is given by 

f(0) = cqI , Co = —8 sinh^^ r] cosh rj sinh «_ sinh a+ cosh (3- cosh j3+ . (2.6) 

For r] 7 ^ ?vr/2, the Hamiltonian (j1 .If) is related to the hrst derivative of the transfer matrix 
at M = 0, 

= ciF(0)+C2l, (2.7) 


where 


Cl = 


C2 = 


— (l6 sinh^^ ^ 7 ] cosh 7 ] sinh sinh cosh P_ cosh 

sinh^ 7] + N cosh^ r] 

2 cosh 7] 


-1 


( 2 . 8 ) 
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For the special case rj = m 12 (i.e., p = 1), 

t'(0) = dol , do = (—sinh a- sinh a+ cosh /?_ cosh f3+ , (2.9) 

and the Hamiltonian (HID is related to the second derivative of the transfer matrix at m = 0 

CH], 

H = dit"{0) , di = (—(32 sinh a_ sinha+ cosh/3_ cosh/3+)~^ . (2-10) 


3 Case I: p even 


For Case I dUD, the Hamiltonian (HID becomes 

H = Hq + - (coshp tanh/3_(T^ + sech — coshp tanh+ sech l3+a^) , (3.1) 

which is Hermitian for (3± real. The eigenvalnes A(m) of the transfer matrix (j2.H] are given 

by [Hj 


A(m) = h{u) 


Q{u + pr]) 
Q{u) 


+ h{-u - rf) 


Q{u- pt]) 
Q{u) 


(3.2) 


where ^ 


h{u) 


A ■ i2V+i/ , xSmh(2M + 2p) 

= 4 smm^'^^^ (u + p) -—-A sinh(M — p) 

sinh(2M + p) 

X (cosh u + i sinh /?_) (cosh u — i sinh /?+) , 


(3.3) 


and 


M 

Q{u) = n sinh 




sinh 


iu + Uj + p'^ , 


M = N + p. (3.4) 


The zeros Uj of Q{u) satisfy the Bethe Ansatz eqnations 


h{uj) _ Q{uj-pp) 

h{-Uj - p) Q{uj + pp) ’ 


j = l,... ,M. 


(3.5) 


^We find that the function h{u) given by Eq. (12) of the Addendum m, to which we now refer as 
hoid{u), leads to p — 1 “Bethe roots” which actually are common to all the eigenvalues, and which therefore 
should be incorporated into a new h{u). In this way, we arrive at the expression which is equal to 

hoid{u) ; and at the M value in 13.411 . which is equal to Mold — {p — 1), where Mou is given by Eq. 

(13) of the Addendum [n|. 
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More explicitly, in terms of the “shifted” Bethe roots Uj = Uj + ^, 


sinh(Mj + 
sinh(-Uj — 


ri\ N 2N+1 


sinh(2'Uj + r]) sinh(Mj — 
sinh(2-Uj - r]) smh.{uj + 


2 ^^ -- y-.^j -ij , 2 

'cosh(Mj — |) + isinh P-\ f cosh(-Uj — |) — isinh/?+ \ 
cosh('Uj + ^) + i sinh 13- ) \cosh(-Uj + |) — i sinh /?+ / 

cosh(i(Mj -Uk + v)) cosh(i(Mj + Uk + v)) . _ ^ 

cosh(i(Mj - Mfc - 7])) COsh(i(Mj + Mfc - 7])) ’ ^ 




(3.6) 


The energy eigenvalues are given by (EZD 


E — ciA^(O) + C2 


= Cih(O) 


Q{pv) 


h'(0) ^ Q'ipv) Q'(0) 


Q(0) [h(0) Q{pr]) Q(0) 


+ C2 . 


Using the fact 




where the second equality follows from (1231), we arrive at the result 


M 


E = — sinh^ ri — 
9 ' 


• w- I ^(iV " 1) COsh 7] . 

^ Sinh(77j sinh(77j + ^) 2 


(3.7) 


(3.8) 


(3.9) 


Numerical investigation of the ground state for small values of N and p (along the lines 
of jS]) suggests making a further shift of the Bethe roots, 


— '^j 


in p in 

Y = “i + 2 “ Y 


(3.10) 


in terms of which the Bethe Ansatz Eqs. (EnD become 


2Ar+2 


cosh( 77 j + 2 ) ^ sinh(77j + |) cosh(77j — ^■) 


cosh (77 


2' 
9 2' 


sinh ( 77 j — |) cosh(77j + 


X 


sinh(|(Y + /3_ - 2)) cosh(|(77j - (3- - f)) sinh(i(77j - /3+ - |)) cosh(|(77j + /3+ - ^)) 
sinh (i (77 j + j3- + ^ 


. . __ .2^7 _ . _ 

;(77j + f3- + f)) cosh(i(77j - /3- + ^)) sinh(|(77j — f3+ + f)) cosh(i(77j + /3+ ^ 


cosh(|(77j -Uk + p)) sinh(|(77j +77^ + p)) 
cosh(|(77j - 77fc - p)) sinh(|(Y + Uk - p)) ’ 


J 
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Moreover, we find that for suitable values of the boundary parameters (3± (which we discuss 
after Eq. (jd.ddjl below), the N + p Bethe roots {hi,... , u^+p) for the ground state have the 
approximate form ^ 


f ± f 

I + ZTT , 




N 
’ 2 
p ’ 
’ 2 


(3.12) 


where {vj , are all real and positive. That is, the ground state is described by y “strings” 
of length 2, and | pairs of strings of length 1. 


We make the “string hypothesis” that ()3.12j) is exactly true in the thermodynamic limit 
{N —>■ cx) with p hxed). The number of strings of length 2 therefore becomes inhnite (there 
is a “sea” of such 2-strings); and the distribution of their centers {vj} is described by a 
density function, which can be computed from the counting function. To this end, we form 
the product of the Bethe Ansatz Eqs. (ITTTD for the sea roots Uj ± y. The result is given by 


[e3(Aj)^5fi+j2b_(Aj)5'i_j2b_(Aj)5'l+i26+(Aj)fi'l-i26+(Aj)] 


-1 


'n/2 

2 

2 p/2 

J^e2(Aj — Afc)e2(Aj -I- A^) 

n n [^2(Aj - xl^'’)g2{Xj -h A^“^) 

k=l 

a=l k=l 


where we have used the notation (Uni), as well as 


Vj pXj , 

and (see HOI and references therein) 
sinh (/i(A + f)) 


= hA 


f'n(A) 


9„(A)=e„(Ad=“) = ^2TkT±M 

^ ^ ^ 2p^ cosh(p(A-f)) 


(3.13) 


i = 1) • • •) IT 


N 

’T’ 


(3.14) 


(3.15) 


sinh (p(A - f)) ’ 

Taking the logarithm of (EH, we obtain the ground-state counting function 
h(A) = ^|(4iV + 4)gi(A) + 2ri(A) - 2qs{X) 

-ri+i2b.{X) - ri_i26_(A) - ri+i26+(A) - ri_i2b+(A) 

N/2 2 p/2 

~2 ^ [g2(A - Afc)g2(A + Afc)] - ^ ^ |^r2(A - A[,“^)r2(AA^“^) |, (3.16) 


k=l 


a=l k=l 


where qn{X) and r„(A) are odd functions dehned by 

Qn{X) = vr-1-/ Inen(A) = 2 tan“^ (cot(n/i/2) tanh(/iA)) , 


r„(A) = i\ngn{X). (3.17) 

"^Due to the periodicity and crossing properties Q{u + 2iTr) = Q{—u — q) = Q{u), the zeros Uj are defined 
up to Uj 1-^ Uj + 2i7r and Uj —Uj — rj, which corresponds to Uj i—> Uj + 2iTT and Uj i—> —Uj — iir, respectively. 
We use these symmetries to restrict the roots to the fundamental region Ke Uj > 0 and —it < 3m Uj < tt. 
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Noting that 


N/2 

[<l2{^ — ^k) + (l2{^ + -^fc)] 

k=l 


N/2 

^ 5'2(A — Afc) — g2(A), 

k=-N/2 


(3.18) 


where X-k = —Xk, and letting N become large, we obtain a linear integral equation for the 
ground-state root density p(A), 


P(A) 


1 dh /■“ 

-- = 2a,{X)-J dX'a^iX - X') p{X') 

^^|4ai(A) -|- 26i(A) — 2a3(A) -I- 2a2(A) — &i+ 2 ife_(A) — &i- 2 ife_ 

2 p/2 

bl+2ib+{X) — &l-2ib+(A) — [^2 (a - A^“^) -h &2 (a -h A^“^)j 


(A) 


(3.19) 


(3.20) 


where we have ignored corrections of higher order in 1 /iV when passing from a sum to an 
integral, and we have introduced the notations 


(\\ - —— M^ _ P sin(n/x) 

” 271 dX^"' TT cosh(2/iA) — cos(n/i) ’ 

h (\) = ^ ^ r (X)- ^ 

” 271 dX ^ TT cosh(2/iA)-|-cos(n/i) 


(3.21) 


Using Fourier transforms, we obtain ® 


P(A) — 2 s(A) -h ]^-^(A) > 


(3.22) 


where 


1 ^OO 1 1 

s(A) = — / doj - , , , , =-—— , 

27r 7-00 2cosh(a;/2) 2cosh(7rA) 


(3.23) 


and 


R{uj) = 


2 (1 - 1 - d2{uj)) 


4 ai(a’) -|- 26 i(a;) — 203(0;) -|- 202(0;) — 6i+2i6_ (t^^) ~ &i-2ib_ (1^) 


2 p /2 


bl+2ib+{^) — bl-2ib+{^) — 2EE cos(o;A^“^) 62 ( 1 ^) \ , 

a=l k=l 


(3.24) 


®Our conventions are 


/(o) = / /(A) dX , /(A) = — / /(o) do. 
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with 


an{uj) = sgn(n) 
bJuj) = 


sinh ((z/ — \n\)ui/2) 
sinh {i'uj/2) 

sinh {nuj/2) 


0 < |n| < 2z/, 
0 < 3?e n < 1 /. 


(3.25) 

(3.26) 


sinh {-uui/2) ’ 

Tnrning now to the expression (El for the energy, and invoking again the string hypoth¬ 


esis (EH, we see that 

. M 


E = — sinh^ 7] 


—^^-1 —(N — 1) coshn 

^ cosh(Mj — 2) cosh(Mj -f 2) 2 


^ V/2 

— sinh^r?!—2^^ — 
2 ^ si 


2 p/2 

EE 


^ smh{vj - f) smh{vj + §) ^ ^ cosh(nJ“^ - f) cosh(n5“^ + 2) 


-|--(A^ — 1) cosh7] 


271 sin fi 


V/2 2 p/2 

{ ^ai(Aj) + 2 

j=l a=l j=l 


COS /i . 


(3.27) 


Repeating the manenver (ITTHD in the snmmation over the centers of the sea roots, and 
letting N become large, we obtain 


E = 


TTSin p 


/i 


TT Sin p 
/i 


N/2 2 p/2 

a=l j=l 


j=-N/2 


2 p/2 


N I dX ai(A) p(A) — ai(0) -|- EE'’i(^f)} + )(^-i)<=“f‘' 

a=l j=l 


where again we ignore corrections that are higher order in 1/N. Snbstituting the resnlt 
(Km for the root density, we obtain 


E Eifiiiii Ef^QiiYidary i 

where the bnlk (order N) energy is given by 

2iV7rsinp 


(3.29) 


^bulk 


/i 


dX ai(A) s(A) -|- -Ncos/i 


= —Nsin^/i / dX 


[cosh(2pA) — cos /i] cosh(7rA) 2 


+ -Ncosfi, (3.30) 














which agrees with the well-known resnlt ng. Moreover, the boundary (order 1) energy is 
given by 


2 p/2 

TT sin U ( T , / \ fa) \ 1 1 

^boundary / j / j ^l\^j ) j COS /i , 

^ n=l 


(3.31) 


where I is the integral 


I = 


d\ ai(A) [RiX)-5{X)] = 


27r 


du ai(a;) R{uj) — 1 


2 p/2 


a=l j=l 


duj s(a;)< 2ai(a;) -|- hi{uj) — d^ioj) — 1 


2 L 


^l+2i6_ (l^) + bi-2ib_ (co’) -|- (co’) + (l^) 


(3.32) 


where we have used the fact s{u!)b 2 {uj) = bi{u). Remarkably, the dependent contri¬ 
bution in (HSH) is exactly canceled by an opposite contribution from the integral I fl3.32jl . 
Writing the boundary energy as the sum of contributions from the left and right bound¬ 
aries, Eboundary = + ^boundary^ ^e couclude that the energy contribution from each 

boundary is given by 

sinp , 1 j'sinh((z/— 2)a;/4) 1 

boundary “ ^ J^ 2 cosh(a;/2) I 2 sinh(z/a;/4) 2 

sinh(a;/2) cosh((i/— 2)a;/2) sinh(a;/2) cos(6±a;)'i 1 (o oo\ 

sinh(z/a;/2) sinh(z/a;/2) / 

One can verify that this result coincides with the result (EHD with a± = 1. As shown in the 
Appendix, the integrals in ()3.3()|1 and ()3.33j) (with p even) can be evaluated analytically. 


We have derived the result for the boundary energy under the assumption that 

the Bethe roots for the ground state have the form which is true only for suitable 

values of the boundary parameters P±. For example, the shaded region in Fig. C]denotes the 
region of parameter space for which the ground-state Bethe roots have the form (I3.12|) for 
p = 4 and N = 2. For parameter values outside the shaded region, one or more of the Bethe 
roots has an imaginary part which is not a multiple of 7i/2 and which evidently depends on 
the parameter values (but in a manner which we have not yet explicitly determined). As p 
increases, the hgure is similar, except that the shaded region moves further away from the 
origin. 


A qualitative explanation of these features can be deduced from a short heuristic argu¬ 
ment. Indeed, let us rewrite the Hamiltonian (ED) as 


n 


Ho + hlal + + h%a 


X 

N ? 


(3.34) 
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Figure 1: Shaded region denotes region of the (/3+ ,/3_) plane for 
which the ground-state Bethe roots have the form (ITT^ 
for p = 4 and N = 2. 

where the boundary magnetic helds are given by 

1 1 
~ 2 ^ fanh /?_ , ~ 2 

= --cosh r/t anh, h^ = -sech/9+. (3.35) 

For ^ 0 (i.e., the shaded regions in Fig. 1), the boundary helds in the x direction 

are small; moreover, hlh^j < 0; i.e., the boundary helds in the z direction have antiparallel 
orientations, which (since N is even) is compatible with a Neel-like (antiferromagnetic) 
alignment of the spins. (See Fig. | 21 ) Hence, the ground state and corresponding Bethe roots 
are “simple”. 


1 I-** fl 




Figure 2: Antiparallel boundary helds (big, red) are compatible 
with antiferromagnetic alignment of spins (small, blue) 

On the other hand, if \P±\ are small (the unshaded region near the origin of Fig. 1), then 
the boundary helds in the x direction are large. Also, if /3+/3_ < 0 (the second and fourth 
quadrants of Fig.l, which are also unshaded), then > 0; i.e., the boundary helds in the 
direction are parallel, which can lead to “frustration”. (See Fig. IHl) For these cases, the 
ground states and corresponding Bethe roots are “complicated”. 
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Figure 3: Parallel boundary fields (big, red) are not compatible 
with antiferromagnetic alignment of spins (small, blue) 

4 Case II: p odd 


For Case II (EZI, the Hamiltonian (dD becomes 


= Tto + - sinh rj (cosech a_ai + cosech o+cr^) . (4.1) 

We restrict a± to be purely imaginary in order for the Hamiltonian to be Hermitian. We 
use the periodicity a± ^ a± + 27ii of the transfer matrix to further restrict a± to the 
fundamental domain —tt < a± < n. The eigenvalues A(m) of the transfer matrix (I2.1j) 
are given by na ' 


where ^ 
h^^\u) 
6{u) 


A(n) 


5{u-vi) Q2 {u-7]) 6{u) Q 2 {u + r]) 

h^^^u — T]) Qi{u) hP)(n) Qi{u) 


— rj) 


Q2{u) 


+ h^'^\u) 


Qiiu + rj) 
Q2 {u) 


1 


(4.2) 


+ 2„) oosli-(« + ,,) oosh(« + 2„) _ ,,W,j , _ ,4,3) 

sinh(2M + 3?]) 

{u)h^‘^^ (u) sinh(M + rj + a_) sinh(n + rj — a_) sinh(M + rj + «+) sinh(M + rj — «+) , 


and 

Ma 

Qa{u) = sinh(M — nj“^) sinh(M + + ??) , a = 1,2 , (4.4) 

j=i 

with 

Mi = ^(7V + p + l), M2 = ^(iV + p-l). (4.5) 

®The function Q 2 {u) here as well as its zeros are shifted by rj with respect to the corresponding 

quantities in to which we now refer as “old”; i.e., Q 2 {u) = Qf‘^{u — rj) and — rj. 

^Similarly to Case I, we hnd that the functions h'^^\u) given in Eqs. (A.5), (A.6) of ^21 Fad to “Bethe 
roots” which actually are common to all the eigenvalues, and which therefore should be incorporated into a 
new h^^\u). In this way, we arrive at the expression for h^^\u) in and the corresponding Ma values in 
( 14 . 511 . 
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As remarked in the Introduction, the expressions for the eigenvalues ( 03 ) correspond to gen¬ 
eralized T — Q relations (O. For generic values of a±, we have not managed to reformulate 
this solution in terms of a single Q{u). The zeros of Qa{u) are given by the Bethe 

Ansatz equations, 


S{uf^) — rf) 

Q 2 {uf -ri) 

2 = 1,2,.. 

• ,Mi, 


5{u^P — rj) h^^'){uf'’) 

Q 2 {uf^ + v) ’ 


1 

1 

Qi{uf + v) 

2 = 1,2,.. 

. ,M2. 

(4.6) 


Qi{uf -7]) ’ 


In terms of the “shifted” Bethe roots the Bethe Ansatz equations are 

J J ^ 


sinh(-u|.^^ -h f) 


2V+1 


cosh(M|-^^ - I) 


sinh(Mj^^ “1)7 cosh(Mj 


sinh(M7^ -f a — |) sinh(-Uy^'’ — a 


X- 


+ f) 

( 1 ) 


3 


sinh(M7^ — a- + ^) sinh(-Uj-^'’ + a- + ^) sinh(Mj-^'' — a+ + ^) sinh(-Uj-^'' -|- a+ 


',(1) 


( 1 ) 


-,(i) 


"3 I 2^ I I 2, 

sinh(-u7^ — + rf) sinh(-u7^ -f- + rj) 


n sinh(F5.^) 


u 


( 2 ) 


rj) sinh(Mj.^^ -|- — rj) 


j = !,■■■ ,Mi. 


(4.7) 


and 


/sinh(Mj^^ + 2 )^ cosh(M® - 2) 

\sinh(M® - 2 ) / cosh(Mj-^^ -f |) 

sinh(Mj-^^ - -f r]) sinh(Mj-^^ r]) 

sinh(-u® - - Tj) sinh(-uj^^ u'^"’ - rj) 


(4.8) 


The energy is given by 


2 Ma 


E = — sinh^ T] EE 


^ sinh(M|“^ - 2) sinh(-u)“^ -1- 2) 


(a) 


-I- -(A^ — 1) cosh? 7 . 


(4.9) 


Indeed, for p > 1, we obtain this result by following steps similar to those leading to (uni). 
For p = 1, we use (O and (ICTl) instead of (ESD and (EZZI; nevertheless, the result (03 
holds also for p = 1. 

From numerical studies for small values of N and p, and for suitable values of the bound¬ 
ary parameters a± (which we discuss after Eq. fj4.2t)|l below), we hnd that the ground state 
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is described by Bethe roots and of the form ® 


,(1,1) 


• 7=12 — 


( 1 , 2 ) ^ . 7 = 12 2±1 
iJj • J ^ ^ y ■ ■ ■ y 2 


,( 2 , 2 ) 


( 2 , 1 ) 


• 7 = 12 — 


J_ IZL • j = I 2 ^ 

'Of • J ^ 1 ^ 1 ' ' ' 1 0 


U ,(4.10) 


respectively, where are all real and positive. 


We make the “string hypothesis” that (I4.1()|l remains true in the thermodynamic limit 
(iV —>• cxo with p hxed). That is, that the Bethe roots for the ground state have the 

form 




^3 

(a,2) 


(a,l) 


: j = 1,2 ,... , M(„4) 


+ Y) : J — 1)2,..., M(a^2) 


a = 1,2, 


(4.11) 


where are all real and positive; also, = M( 2 ,i) = ^y and M(i^ 2 ) = ^^y 1 ^( 2 , 2 ) = 

2^. Evidently there are two “seas” of real roots, namely and 

We now proceed to compute the boundary energy, using notations similar to those in 
Case I. Defining 


the Bethe Ansatz equations (j4.7|l . (14.811 for the sea roots are 


(A;.‘’‘')e,-2«_ (A'‘'pe,_2„, (A'*’*') 


,(14)) 


,( 1 , 1 ) 


(1,1)a 


,(1,1)a 


,( 1 , 1 )) 


n -1 


(4.12) 


(4.13) 


N/2 




( 2 , 1 )) 
k ) 


(p-l)/2 


n k2(Af’^)-Af^)),2(Ar^+Ar) 


,(1,1) 


,( 2 , 2 )) 


and 


k=l 


e.(Af‘>)“+‘ji(A, 


k=l 


( 2 , 1))-1 


N/2 (p+l)/2 

n h(Af- A<‘A>)e2(Af ■'+ A(‘A))1 Yl U(Af‘> 


(4.14) 

A7b92(Af A) + aAA')1 . 


k=l 


k=l 


respectively, where j = 1,... , y- The corresponding ground-state counting functions are 
(A) = — I {2N + l)qi (A) — ri (A) — gi+ 2 a_ (A) — gi- 2 a_ (A) — gi+ 2 a+ (A) — q'i- 2 a+ (A) 

N/2 (p-l)/2 

“ ^ ['?2(A - A^^’^^)-h g2(A 4-A[,^’^^) - ^ |^r2(A - A[,^’^^)-1-r2(A-4 A^^’^^) |(4-15) 

k=l 


k=l 


,(° 


The periodicity and crossing properties of Qa{u) imply that the zeros are defined up to 

+ ITT and 1 —*■ respectively. We 


+ iiT and 77^“^ 


— rj, which corresponds to 12^°^ 


use these symmetries to restrict the roots to the fundamental region (Re > 0 and — § < ^rn 
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and 


h'"nA) = ^{(2/V + l)9i(A)-r,(A) 


“ ^ [92(A - ^ ^^(A + A[^a)j 


(P+l)/2 

E [r2(A - A<‘’^') + r,(A + A<‘’^>) 



k=l k=l 

Repeating the manenver (ITTHl) in the snmmations over the sea roots, and letting N become 
large, we obtain a pair of conpled linear integral eqnations for the gronnd-state root densities 
p(“)(A), 


1 dh(l) 

'iA'a2(A-A')pW(A') 


+ ®2(A) — bi{X) — ai+2a_ (A) — ai-2a_ (A) 

(p-l)/2 


fli+2a+(A) — ai_2a+(A) — |^&2(A ~ A^ ’ ^) + 62(A + A^ ’ ^) I, ( 4 . 17 ) 


k=l 


and 


p(^)(A) = 


1 

N~dr 

1 r 


= 2 ai(A) - / dX' a 2 {X - A') p^^\X') 


(P+l)/2 


k=l 


+ + «2(A) - 6i(A) - ^ |^62(A - + 62(A + A[,^’^^) | .( 4 . 18 ) 


It is straightforward to solve by Fonrier transforms for the individnal root densities. However, 
we shall see that the energy depends only on the snm of the root densities, which is given by 


pW(A)+p'"'(A)= 4 s(A) + ^K(A), 




where 


n{uj) = = 


2 ai(a;) + 202(0;) — 26 i(a;) — ai+2a_(o;) — hi-2a_(o;) 


( 4 . 19 ) 


1 + 02 ( 0 ;) 


2 ^(<2,2) 


- ai+2a+(o;) - di-2a+{uj) - 2^ ^ cos(a;A^“’^^) 62(0;) 


( 4 . 20 ) 


a=l k=l 


The expression (Oil for the energy and the string hypothesis (OUl) imply 

2 N/2 


E = --sinh^T^I 


Tt ^ sinh(n^“’^^ - 2) sinh(n)“’^'' + 2) 


(a,l) 


a=i j 
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2 ^{a,2) 


a'^i j=i cosh(vj“’^^ - 2) cosh(t;]“’^^ + f) 

2 N/2 2 ^^(0,2) 


+ -(A^ — 1) coshr; 


TTSin /i 

/i 


{eE“'(^5“’“’)+e e ('.(Ar’jj+jiiv-i) 


cos /i 


a=l j=l 
2 Af/2 


a=l j=l 


2 ^(a,2) 


^{iE E a.(Ar')-a.(0)+5:E'>.(r>)}+>-l) 

' n — 1 A \T/0 n — 1 o — 1 


l2^ ^ 

a=l j=-Nl2 a=l j=l 

r ‘'A «.(A) [p<'>(A) + pl='(A)] - «.(0) 

^ ^ J —oo 


cos /i 


2 ^(a,2) 

E E '>.(A<“'^’)} + 5(iv-i) 


cos /i . 


(4.21) 


a=l i=l 

Substituting the result (ICTD for the sum of the root densities, we obtain 

^ Ebulk 4“ ^boundary i 


(4.22) 


where the bulk (order N) energy is again given by (jd.dOjl . and the boundary (order 1) energy 
is given by 


^boundary 

where X is the integral 

^ roo 

I = - i 


TTSin /i 
h 


2 *^(a.2) 


{i+EEMAr>)}-i 


cos Id . 


(4.23) 


a=l j=l 


dX ai(A) [7^(A) - 2(5(A)] = 


47r 


duj hi (a;) TZ^u) — 2 


2 ^(a,2) 

EEMAr’)+^/ 

_ 1 A 1 d — c 


duj s(a;)< ai(a;) — 6i(a;) — 1 


a=l j=l 


'2 [®l+2a-(‘^) + ®l-2a_(<^) + hi+2a+(<^) + hi_2a+(<^)] 


(4.24) 


Once again there is a remarkable cancellation among terms involving Bethe roots which are 
not parts of the seas, namely, Writing Eboundary as the sum of contributions from the 

left and right boundaries, we conclude that for the parameter values ® 


1,1 1 

- < a± < Id -, 

2 - I ±1 2 ’ 


(4.25) 


^This restriction arises when using the Fourier transform result to evaluate ai+ 2 a±(w) + ai_ 2 ai-(w). 
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the energy contribution from each boundary is given by 


sin/i , 1 [ cosh((i/— 2)a;/4) 1 

boundary = ^ J 2 COsh(o;/2) I 2 COsh(z/w/4) 2 

+ «toMW2)cosh((,.-2|a^|V/2)|_l^^ 

sinh(i/a;/2) J 4 

This result agrees with the result m with b± = 0 and a± values (|03)- As shown in the 
Appendix, the integrals (with p odd) can also be evaluated analytically. 

We have derived the above result for the boundary energy under the assumption that the 
Bethe roots for the ground state have the form ()4.1()|1 . which is true only for suitable values 
of the boundary parameters a±, namely, 

< |a±| <, a+a_>0, (4.27) 

where u = p + 1. For parameter values outside the region (103), one or more of the Bethe 
roots has an imaginary part which is not a multiple of 7r/2 and which evidently depends on 
the parameter values. One can verify that the region (R37I1 is contained in the region (jOSl)- 

As in Case I, it is possible to give a qualitative explanation of the restriction ()4.27j) by a 
heuristic argument. Indeed, let us rewrite the Hamiltonian dO} as 

n = no + , (4.28) 


where the boundary magnetic helds are given by 


h 


X 

1 


- sinh p cosech a_ , 


Hn — 


- sinh p cosech . 


(4.29) 


For a± ~ 171/2 or —177/2 (i.e., a± ~ v/2 or —v/2\ the boundary magnetic helds in the x 
direction are small and parallel. Hence, the ground state and corresponding Bethe roots are 
“simple”. Outside of this region of parameter space, the boundary helds in the x direction 
are large and/or antiparallel, and so the ground state and corresponding Bethe roots are 
“complicated”. 


5 Discussion 

We have investigated the ground state of the open XXZ spin chain with nondiagonal bound¬ 
ary terms which are parametrized by pairs of boundary parameters, in the thermodynamic 
limit, using the new exact solutions ocn] and the string hypothesis. This investigation 
has revealed some surprises. Indeed, for Case I dini), the ground state is described in part 
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by a sea of strings of length 2 (IXT^ . which is characteristic of spin-1 chains [20] • For Case 
II (ini), the energy depends on two sets of Bethe roots (Q, and in fact on the sum of the 
corresponding root densities dOII). For each case, there is a remarkable cancellation of the 
energy contributions from non-sea Bethe roots. 

Perhaps the biggest surprise is that, for the two cases studied here, the boundary energies 
coincide with the result (EHD for the constrained case (d, even when that constraint is not 
satished. This suggests that the result m may hold for general values of the boundary 
parameters. A hrst step toward checking this conjecture would be to extend our analysis to 
the “unshaded” regions of parameter space, where the ground state has Bethe roots whose 
imaginary parts depend on the parameters. 

We have not computed the Casimir (order 1/A^) energy for the two cases ()l.f)|l . (II. 7j) . The 
computation should be particularly challenging for the former case, due to the presence of a 
complex sea. It would be interesting to investigate excited states, and also applications to 
other problems, including the boundary sine-Gordon model. We hope to be able to address 
such questions in the future. 
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A Appendix 

Here we present more explicit expressions for the bulk and boundary energies. 

A.l Case I: p even 

The integral appearing in the bulk energy (HOnil for p even is given by 


/i 



E 



(A.l) 
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The parameter-dependent integral appearing in the boundary energy is given by 


T (r) = ^ r Hi ; sinh(a;/2) cos(3:a;) 

^ 27r 2 sinh(z/a;/2) cosh(a;/2) 


i=i 


where the function 6„(A) is dehned in fjd.21|l . Moreover, 


1 , sinh(cu/2)cosh((z/-2V/2) ^ 

^ J_^ 2sinh(i/o;/2) cosh(o;/2) “ ’ 


and 


1 

271 


sinh((z/-2V/4) _ 

4sinh(i/a;/4) cosh(c<;/2) 2 ^ ^ 


It follows that the boundary energy 


is given by 


where Ji and hix) are given by (EH) and (Ha, respectively. 


A.2 Case II: p odd 


(A.2) 


(A.3) 


(A.4) 


Etundary = + hi^P " l)/2) + l2{b±) - ^ ^ COS /i , (A.5) 


The integral appearing in the bulk energy (HOnil for p odd is given by 

pco poc 


h = 


dX ai(A) s(A) = 


27r 


doj hi (a;) s(a;) 




TT^ L 


P-1 

2 


1 + 2/i JCOt(j/i) 
i=i 


(A.6) 


The parameter-dependent integral appearing in the boundary energy (OHll is given by 

j (r\ = — r di ; sinh(o;/2) cosh(xo;) 

^ 27r 2 sinh(z/a;/2) cosh(a;/2) 


Moreover, 


TTsinfxTr') 


p-i 

2 


a: -|- ^^(—1)-^ cot(j/i) sin(2xj/i) 
i=i 


1 /•“ 
27r 


, , cosh((t/ - 2)o;/4) _ 1 

4cosh(i/a;/4) cosh(a;/2) 2 ^ ^ 

18 


(A.7) 


(A.8) 

















It follows that the boundary energy is given by 


+ \m + mp+l - 2|a±|)/2) - i} - icos^, (A.9) 

where Ji and l 2 {x) are given by jH) and dm), respectively. 
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